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Abstract 



We study the properties of the lower bound on the exchange-correlation energy in two dimensions. First we review the derivation 
of the bound and show how it can be written in a simple density-functional form. This form allows an explicit determination of the 
prefactor of the bound and testing its tightness. Next we focus on finite two-dimensional systems and examine how their distance 
i__ from the bound depends on the system geometry. The results for the high-density limit suggest that a finite system that comes as 
^—t close as possible to the ultimate bound on the exchange-correlation energy has circular geometry and a weak confining potential 
IT) with a negative curvature. 
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1. Introduction 
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The lower bound on the quantum mechanical part of the 
Coulomb interaction energy, commonly known as the Lieb- 
Oxford (LO) bound [ 1], is a key concept in many-body physics. 
The bound is not only of fundamental importance for, e.g., an- 
alyzing the stability of matter yfl, but it has also been used 
as a key constraint in the construction of several exchange- 
correlation functionals within density-functional theory (DFT), 
which is a standard tool in electronic-structure calculations of 
atoms, molecules, solids, and various nanoscale systems j^]. 
Recently, substantial efforts have been directed at testing the 
tightness of the original three-dimensional LO bound 14j,|5[] and 
in finding tighter forms [3 |7|] . 

In view of the growing interest in two-dimensional (2D) sys- 
tems, the existence and properties of the 2D form of the lower 
bound on the exchange-correlation energy are of immediate rel- 
evance. In this paper we review the known form of the 2D 
bound and show how it can be derived from scaling relations. 
We discuss the tight prefactor for the bound obtained from the 
properties of the homogeneous 2D electron gas (2DEG) in the 
low-density limit. We test this bound for finite 2D systems, 
and, in particular, by varying the shape of the system for a sim- 
ple two-electron system, we propose a set of general properties 
for a finite 2D system which is as close as possible to the lower 
bound on the exchange-correlation energy. 
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2. Two-dimensional bound 

Lieb, Solojev, and Yngvason |@] (LSY) have rigorously de- 
rived a 2D form of the LO bound, which can be expressed in 
terms of the indirect part of the interaction energy: 



W„m = {v\V„W} - U[n] > -C J d 2 rn 3/2 (r), 



(1) 



where V e . 



■">] i 



is the Coulombic e-e interac- 



tion operator, ^(ri, r#) is any normalized 2D many-body 
wave function, n(r) is the corresponding density, and U[n] is 
the classical Hartree energy. For the prefactor LSY estimated 

C<Clsy = 192V27T*481. 

Interestingly, the exponent 3/2 in Eq. ([TJ follows directly 
from universal scaling properties of the e-e interaction J3]. Un- 
der homogeneous 2D coordinate scaling, r — > yr (0 < y < 
oo) fgj] the (27Y)-dimensional many-body wavefunction scales 
as l^ri . . . r N ) -> ^(Ti . . . r N ) = y 2N/2x ¥(yri . . . yr N ), preserv- 
ing normalization. This produces the number-conserving scaled 
density n(r) -> n r (r) = y 2 n(yr). Further, W^PF] -> W^PF?] = 
yW xc [*F], since both the Coulomb interaction and its Hartree ap- 
proximation scale linearly. Denoting the exponent in Eq. ([TJ by 
X, we then find the relation 



yw xc m>-c r 



d 2 rn x (r), 



(2) 



so that consistency with Eq. ([JJi gives x = 3/2. 

It was also conjectured in Ref. [7] by the present authors that 
the prefactor can be decreased to C = 1 .96 which corresponds 
to a significant tightening of the original 2D bound with Clsy ~ 
481. The tightest bound, i.e., the lowest exchange-correlation 
energy corresponds to the 2DEG in the low-density limit, and 
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hence all other 2D systems (including all real, finite systems) 
are energetically above this bound. However, it remains to be 
examined which type of a finite system is closest to the bound. 

3. Density-functional form of the bound 

Here we reformulate the bound in Eq. ([]]) in terms of den- 
sity functionals. The right-hand side can be written in terms 
of the expression of the local-density approximation (LDA) for 
the electronic exchange in 2D, 

£LDA [n] = -A J d 2 rn V2 (r), (3) 

where A = 4 V2/(3 y/jr). This formula has the same scaling 
with respect to the density as Eq. ([TJ. Note that Eq. is ex- 
act for the exchange energy of the 2DEG (constant n) by con- 
struction [ 10]. The left hand side of Eq. (0 can be written in 
terms of the exchange-correlation energy as defined in DFT, 
i.e., E xc [n] = W xc [n] + T c [n] > W xc [n], where T c is the dif- 
ference between the many-body kinetic energy and the (single- 
particle) Kohn-Sham kinetic energy and it is always positive. 
Now, Eq. (0 becomes 

E«[n] > | E\ Dk [n]. (4) 

For any 2D system, we can now consider the density functional 

E xc [n] C 
A W = -^h^ * 7 (5) 

E^ DA [n] A 

with E xc [n] = E x [n] + E c [n]. As mentioned above, the tightest 
2D bound corresponds to the 2DEG in the low-density limit 
which yields the maximum value of Eq. (0, i.e., /12degI/.s — > 
oo] = 1 .84, where r s = 1/ yjmi is the density parameter in 2D. 

4. Testing the bound 

4.1. General remarks 

In practice, calculation of A[n] for finite systems and thus 
testing the bound is seriously limited by the lack of reference 
data for exact exchange-correlation energies and exact densi- 
ties. An exception is the 2D Hooke's atom, i.e., a parabolic 
(harmonic) quantum dot (QD) with two electrons (N = 2), 
for which analytic solutions are known 111 ill , and which yields 
A[n] « 1.55 as the maximum value Q]. 

It is also possible to approximate A from the exact refer- 
ence data solely for the total energy E tot . This requires exact- 
exchange (EXX) or Hartree-Fock calculations to obtain the ex- 
change energy, which, for closed-shell systems, can be written 
as 

E x [n] = -JdrJdr — , (6) 

where t/f, are Kohn-Sham (or Hartree-Fock) orbitals. The cor- 
relation energy is then obtained as E c = E tol - Ef xx . It should 
be noted that Eq. (0 becomes now an approximation due to us- 
ing both exact densities and EXX densities as the input instead 



of consistently using only the exact densities. Nevertheless, 
tests for few-electron parabolic and square-shaped QDs based 
on this strategy have led to values in the range 1.1 < A[ri] < 1.5 
(Ref. [7]). The results for QDs suggest that the largest value for 
A is obtained with N — 2. 

4.2. Effects of geometry 

To consistently analyze the dependence of A[n] on the system 
geometry, we focus in the following solely on the limit r s — > 

0. This corresponds to the noninteracting situation, since the 
kinetic energy scales as rj 2 and the interaction energy as rj l 
(cf. the opposite limit where interactions dominate and lead 
to Wigner crystallization B12ll ). Now, the correlation energy is 
zero, but Eq. ((5) is still a well-defined quantity having a form 

jr m E x [n] 

Furthermore, we set N = 2 so that the exact exchange energy 
in Eq. (0 can be calculated as a simple integral over the den- 
sity; in fact it is exactly minus half of the Hartree energy. To 
find the density we solve the Schrodinger equation for a non- 
interacting singlet state in the presence of an external confin- 
ing potential V(x, y). In the numerical procedure we use the 
octopus code II 1 311 . Note that upon the condition that the po- 
tential scales homogenously with respect to the scaling param- 
eter y (see above), we can choose any prefactor in the potential, 

1. e., any £ in £V(x,y), in order to mimic the true interacting cal- 
culation in the limit r s — > corresponding to n — > oo. This 
can be seen by considering V(r) in the radial Schrodinger equa- 
tion for N - 2, where iffi(r) oc yjnjr) (i = 1,2). The form 
of the equation under uniform coordinate scaling r — > yr (see 
above) shows that V(r) — > y 2 V(yr) = y 2+a V(r), assuming that 
V(r) scales homogenously with y. Hence, a value for y can be 
found for any prefactor if . In other words, changing the pref- 
actor is equivalent to scaling of the density. And in fact, this 
is the case in all potentials given below, as all of them scale 
homogeneously with respect to y. Moreover, Eq. (0 is now 
independent of y. 

First we consider a circular QD defined by a confining po- 
tential of the form 

VcircularM = |lf, (8) 

and a square-shaped QD defined by 

V square U,y) = |*r + |yr (9) 

Here the parameter a determines the "steepness" of the poten- 
tial. 

Figure shows A in the r s — > limit as a function of a. 
Obviously, the two potentials are the same when a — 2, when 
they actually correspond to the 2D Hooke's atom (dotted line) 
with A as 1.107. At a > 2 we find monotonous decrease of A, 
and the limit a — > oo, which corresponds to the hard-wall case, 
leads to A » 1.0963 and 1.0949 in circular and square QDs, 
respectively. Overall, the circular potential gives higher values 
of A than the square one. 

Interestingly, the highest values for A are obtained at a < 2. 
In this regime, a = 1 corresponds to a cone (pyramid) in a 
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Figure 1 : (color online) Values for A in two-electron quantum dots with varying 
external confining potential in the noninteracting limit. 



Figure 2: (color online) Values for A in elliptic quantum dots in the noninter- 
acting limit. The shaded figures represent schematically the symmetry of the 
confining potentials. 



circular (square) QD. The numerical accuracy of our eigenvalue 
solver limits the investigation to a = 0.5, which actually yields 
the largest A as seen in Fig. [T] Therefore, we assume that the 
maximum value for A could be found in a circular QD at a < 
0.5, when the curvature is negative, i.e., d 2 V(r)/(dr 2 ) < 0. 

In Fig. [2] we consider the r s — > limit in an elliptic confine- 
ment, 



VellipticO,.y) =y6x" 



y 2 //3, 



(10) 



where /3 is related to the eccentricity, although this is not a 
formal definition. Similar confinement has been used in the 
QD studies in Refs. |[l4ll and |15J. We find that A decreases 
as a function of and the highest value can be obtained in a 
parabolic system at /3 — 1 (Hooke's atom). 

5. Conclusions and outlook 

Summarizing the results reported in the previous section for 
finite systems, we have found that in order to find a maximum 
value of A in the r s — > limit, the confinement potential should 
be (i) circular instead of square or elliptic, and (ii) it should 
have a negative curvature. Assuming that these geometric ef- 
fects carry over to r s > 0, and in particular to values for r s 
maximizing A(r s ), and combining these findings with previous 
results for harmonic systems with N > 1 and r s > (Ref. |7[]), 
we advance the following conjecture: The maximum value of 
A infinite 2D systems, and hence the closest possible value to 
the ultimate bound on the exchange-correlation energy can be 
obtained in a two-particle system having circular symmetry and 
a weak confining potential with a negative curvature. 

Similar to what was previously observed in 3D S |(J, 
smaller particle numbers in 2D produce larger values of the 
functional A[n] (Ref. |7j]). This should be compared with the 
behavior of the function A(N), which for any N produces an up- 
per limit on A[n] for all densities integrating to this N. In 3D, 
A(N) is known rigorously to be monotonically increasing 



with N. We expect this to be true also in 2D, but have not proved 
it. In any case, the fact that the upper limit A(N) increases with 
N, while the actual value A[n] decreases is not a contradiction. 
It simply means that the LO bound becomes more and more 
generous as N increases, and, conversely, tighter and tighter as 
N decreases. 

Overall, we conclude from the investigation of the present 
paper that, qualitatively and even semi-quantitatively, 2D sys- 
tems behave similarly to 3D systems with respect to the appro- 
priate LO bound. 

Furthermore, in view of the results for the 2DEG and 2D 
Hooke's atom ||7|] we may assume that a finite system with 
the largest A would have rather uniform density (except at the 
boundaries), and A would be very close to the bound value 
Mdeg = 1-84. Although we cannot rigorously prove these as- 
sumptions for the potential or for the density at large r s , we 
hope that these results will encourage similar geometric stud- 
ies on finite systems where correlations are incorporated. This 
would require using, e.g., the density-functional formalism for 
strictly correlated electrons lfl6ll . or an inversion scheme 1 171. 
to reconstruct the exact external potential and the corresponding 
many-body problem for a given density. 
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